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We study the transition temperature Tc{H) of a quasi-one-dimensional (QID) superconduc- 
tivity which is derived from the quantum effect of an electron motion in a strong magnetic field. 
We calculate Tc{H) of both isotropic and anisotropic superconductivity by taking account of the 
optimal momentum of the Cooper pairs and the effect of higher harmonic terms along second 
conducting axis in the tight-binding model. We find that, although Tc{H) of the spin-singlet 
superconductivity is suppressed strongly by the Zeeman effect, the suppression of Tc{H) is not 
very severe if we take the optimal pair-momentum and the higher harmonic terms into account. 
The obtained Tc{H) for the spin-singlet superconductivity is consistent with the experimental 
results in TMTSF safrs. 
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It is known that the spin-singlet superconductivity is 
destroyed by both the orbital frustration and the Pauli 
paramagnetic effect. Recently, the reentrance of the su- 
perconductivity which is caused by quantum effect of or- 
bital motions along the open Fermi sru;facjeiiiji_a,strong 
magnetic field has attracted interest.ErBB'Q'El'&LP The 
anomaly of the reaistivity in a strong magnetic field has 
been observedB'BEj^ in quasi-one-dimensional (QID) or- 
ganic superconductors (TMTSF)2X (where anion X is 
CIO4 or PFg) when the magnetic field is applied along 
the second conducting axis (6 axis). This is thought to 
be a signal of the superconductivity in a strong magnetic 
field, since the critical temperature Tc{H) exceeds Jiath 
the upper critical field Hc2 derived in GL thecffjtiO 



and Pauh paramagnetic limit Hp[T]= I.84rc[Kl 
In the QID system, it has been knowru' 
that the spin-singlet superconductivity is not destroyed 
completely due to the Zeeman effect by constructing 
the Larkin-0|Mcliinnikov-Fulde-Ferrell (LOFF) supercon- 
ducting statetSllB) in which Cooper pair is formed by the 
electrons (k,^) and (— fe + q, |). The electron (— fe + q, j) 
can be on the down-spin Fermi surface for any (fc, "f ) 
on the up-spin Fermi surface in a ID system by choos- 
ing the appropriate q, which is similar to the nesting of 
the Fermi surface in the spin-density-wave (SDW) case. 
Even in the QID system, the "nesting" condition for 
the LOFF state was iJbnuehlJ^a become perfect in the 
strong iTLamicticLfixJd.BBd'Q These theoretical cal- 
culationsEl'B'BQlj'BlZP for Tc(H) based on the approxima- 
tion that the Fermi velocity is independent of the posi- 
tion of the Fermi aiirface. 

Recently, LebedEO^ has shown that the nonlinearity ef- 
fect of the energy dispersion along the a axis on the Zee- 
man splitting causes the finite upper critical field in QID 
systems. He obtained that the critical magnetic field 
for the LOFF state is ff^OFF _ Q,Q^JtjTbH-p, where ta 
and % are the hopping matrix elements along a and b 



axes, respectively. Applying this result to (TMTSF)2X, 
jjLOFF ^ 4Xesla is obtained, which is smaller than the 
experimentally observed value by Lee et al. (at least 
7T).la'Ll^ This result may suggest that the superconduc- 
tivity in this system is a spin-triplet state which is not 
affected by the Zeeman effect. However, Tc{H) mea- 
sured by Lee et al. does not reveal the reentrant behav- 
ior ^gj^gtgd in the case of spin-triplet superconductiv- 

In this paper, we study Tc{H) of a QID spin-singlet 
superconductor by taking the optimal pair momentum 
q and introducing the higher harmonic terms along the 
second conducting axis in the tight-binding model, which 
has been discussed, Jii the field- induced-spin-density- 
wave rmpW)Ea'B0 and the quantum hole effect 
(QHE)Bi3) in FISDW. We consider both isotropic and 
anisotropic pairing states for spin-singlet superconduc- 
tivity, p*|Ce the symmetry of the pairing is still contro- 
versiaLSEJ) 

We consider the anisotropic tight-binding model in- 
cluding the effect of the Zeeman splitting (we take Ti = 
= Co = I, where co is the velocity of light): 

^ = —2taCos{akx) ~ 2tbC0s{bky) — 2t2b cos{2bky) 

—2tj,b cos{3bky) — 2i4b cos(46fcj,) — 2tc cos{ckz) 

-(t^ibH - /i, (1) 

where ta ^ h ^ \t2b\ > \t3b\ > \t4b\ ^ tc, ct^bH is the 
Zeeman energy for t (J,) spin (tr = +(— )) and /i is the 
chemical potential to give the quarter filled electrons. 

The orbital effect of the magnetic field is treated by the 
Peierls substitution, i.e. Ho = E{k —iV — eA). We 
take the vector potential A as A = {0,0,— Hx). Since 
ta 3> tb, we can linearize the energy dispersion along kx 
for each ky . Then the eigenvalues are given by 



j{akx 



(2) 
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where a — sgn(fcx) refers to the right/left sheet of the 
Fermi surface. Although we linearize the dispersion, we 
consider the ky and cr dependence of the Fermi wave 
number ^ and the Fermi velocity along the a axis, 



2taasin{ak^ ^). Note that A:^ 



and 



do 



not depend on k^ when the magnetic field is applied along 
the 6 axis. For given ky and cr, k\^ ^ is obtained by, 

1 



[{2tiy cos{hky) + 2t2b cos{2bky) 



+ 2t^i, cos{3bky) + 2t4f, cos(46fcy) 
The corresponding eigenstates are 



jk r sr^ 



(3) 

(4) 



where G — eHc, Jn is Bessel function and ■q'^ ^ = 
~2at,/vl^^G. 

If we expand ^ to the second order in /i^ around 
kx ~ i/cF in eq. (1), we obtain 

Eky,a,a ~ vpiakx - kp) - 2tbCos{bky) - 2t2b cos{2bky) 

— 21^1, cos{3bky) — 2t4h cos(46fc.y) — f3tb cos^ (bky) 

+ a^BH{l — (3cos{bky)) — 2tc cos{ckz) — (5) 

where kp and vf = 2iaSin(a/cF) are the Fermi velocity 
and the Fermi wave number in the ID case = 0), 
respectively and (5 = ^/2tb/ta. Lebed has discussed that 
the term proportional to j3 in the Zeeman energy causes 
the finite critical field in QID systems. In the following, 
we do not use the expansion of ^ around i/cp (eq. 
(5)). We calculate Tc(ff) by using eqs. (2)~(4). 

The one-particle Green's function in the mixed repre- 
sentation is 



G'^{x,x',ky,k^-UJ,,) = ^e''^(ck^-Gx)+in'(ck^-Gx') 
n.n' 

where G'^{kxTky,uJn) = l/(iw„-e^^ ^.^ anda;„ = (2n+ 
l)7rT is a Matsubara frequency. 

We first study Tc{H) of the isotropic superconductiv- 
ity caused by the on-site attractive interaction A in the 
mean field approximation. The linearized gap cjquation 
for an isotropic superconductivity is obtained ag 



A(.) = AT E E 



dx'A{x') 



x — x' I >(/ 

X ^ ky ^ kz] iOn)G^{x^ X , ky^ 



-^n), (7) 



where d is the cutoff. 

The solutions of the gap equation (7) are written as 



AQ(x)=e"^-}_^A^,e'^''-^ (8) 
I 

where Bloch wave vector Q is taken as —G < Q < G. 



Then eq. (7) is written as a matrix equation 

^^i-^T.<2i'^%^ (9) 

where 

^21,21' = T.T.^U'ih)KkyiQ + NG). (10) 

N ky 

The coefficients S[^i,{ky) for an isotropic superconductiv- 
ity are defined by 

n 

X Jn-l+N(r]Z^^s)Jn-V+N(Tll^^g)- (H) 

In the above Kky{qx) is given by 

Kky{qx) = E G'^{kx,ky]UJn)Gg{qx - kx,-ky;-LJn) 



E 



F I F 

ky .(7 ky ,a 



2 4i7rr' 



(12) 



where 



2vl 



7/K„.a + '"L.a)^ 7 is the expo- 



^k s ' which is the wave number of 



nential of the Euler constant, is the cutoff energy and 

4' is the digamma function. If fc^^ ^ — fc^^ ^ -|- aqx = is 

satisfied, Kky{qx) diverges logarithmically as T goes to 

zero. This logarithmic divergence survives over the ky 

summation in eq. (10) only when Fermi surfaces for the 

up and down spins are "nested". This is not the case 

when the ky dependence of the Fermi wave number is 

taken into account. The "nesting" of the Fermi surface 

becomes worse as H increases, and as a result the critical 

magnetic field has a finite value 
Ineq. (12),fcF 

the Cooper pair of electrons on the left Fermi surface of 
down-spin and right Fermi surface of up-spin, gives the 
information of the "nesting" condition as shown in Fig. 
1(a). InFig. 1(b), we plot a(fcf^_^-fcF_^) for Wia = 0.1, 
tib/tb = and some values of higher harmonic terms 
t2b and tab at H = 5T. We write the x-component of 
"nesting vector" or the wave number of the Cooper pair 
q as qx{s) = 2sij,bE[/vf, and plot aqx{s) in Fig. 1. 

In Fig. 2, we plot Tc{H) of the isotropic supercon- 
ductivity obtained from eq. (9). In the following, we 
take parameters as 2ta = 1950K, tb/ta = 0.1, t4b/tb = 0, 
Tc(0) = 1.35K. The maximum value of Tc{H) for each 
H is obtained by optimizing qxi-s)- 

The lines with squares, solid circles and open circles in 
Fig. 2 are obtained for tc/ta = 0, which corresponds to 
no orbital effect. The optimized Tc{H) in the absence 
of higher harmonic terms is plotted by the line with 
solid circles. We find that the critical field at T = 0, 
is more than two times larger than the 
PauH paramagnetic limit Hp = 1.84Tc(0) ^ 2.5T. This 
result is also larger than that of the squares, where the 
nesting vectoii.q is fixed to be 2^bH/vf, (s — 1) as stud- 
ied by Lebedla) {H^^^^ = O.^y^tJTbHp ~ 4.7T). By 
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Fig. 1. (a) Schematic QID Fermi surface in kx-ky plane in a mag- 
netic field, (b) a{k^ a ~ „) ^ ^ function of ky/{b-K). We 
take parameters as %/ta = 0.1, t46/ti, = 0, quarter filled band 
and H = STesla. The long-dashed (dotted) line represents the 
x-component aqx{s) = 2.safi^H /vp of the "nesting" vector aq in 
s = 1.0(0.9). 



adding the higher harmonic terms t2b/tb = —0.1 and 
tsb/tb — —0.07, i/p*-"^^ becomes larger (the open cir- 
cle line in Fig. 2, 9T). The enhancement of 
Tc{H) and Hp'^^^ due to the optimization of q and the 
higher harmonic terms (t26 < 0) can be understood by 
the "nesting" of the Fermi surface as follows. At low tem- 



-aqx\ < 



peratures, only electrons with vf ^\k^ a~^k 
(5 ~ T can contribute to forming Cooper pairs. The area 
on the Fermi surface, where v^^ Jkl^^^-kl^^^ + aq^l < S, 
is proportional to VS as 5 ^ if down Fermi surface 
touches the up Fermi surface by the translation of q, 
while it is proportional to 5 if two Fermi surfaces cross 
by the translation. Since > (5 for 5 -> 0, H^O'^^ 
becomes largest when the Fermi surfaces touch by the 
translation, which is equivalent to the case when qx{s) 
touches ^ — kf g. in Fig. 1(b). The enhancement 
of Hp'~'^^ ^P-?-^ ^^^^ mechanism has been studied by 
Shimahara.ES' These are two possibilities of choosing op- 



qx{s) touches at ky « 0. If we take the positive t2b, the 
"nesting" becomes worse at fcy « 0, although it becomes 
better at ky « zLir/b as shown by the dot-dashed line 
Fig. 1(b). In this case 

^LOFF 

is not enhanced. We get 



jyLOFF 



the other hand k^ 



5.8T for t2b/tb = 0.1 and t-iblh 

F 



= -0.07. 



k\ g becomes flatter at ky 



On 

« 



when t2b is negative as shown by the dashed line in Fig. 
1(b), resulting in the better "nesting" and larger Hp'^^^ , 
as shown by the line with open circles in Fig. 2. 

Next, we take the effect of orbital motions into ac- 
count. Since the field dependence of the initial slope 
di?^2/dif'|k'c(o) is approximately given by the parame- 
ter tcia"EP in the week field limit, we take parameters 
as 2ta 1950K, 2tc ^ 3.0K, a - 7A, c ~ ISA and 
Tc(0) ~ 1.35K irynrder to fit the initial slope observed in 
(TMTSF)2C104.By) The transition temperature is plot- 
ted as the thick solid line in Fig. 2. This curve seems to 
be consistent with the experiments ip organic supercon- 
ductors (TMTSF)2X by Lee et a/Ml 
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Fig. 2. Transition temperature Tc{H) of an isotropic pairing as a 
function of the magnetic field in the case of tf^/ta =0.1, t4i,/ii, = 
0, Tc(0) = 1.35K for quarter filled electrons. The lines with 
squares and open circles and thick solid line are obtained from 
the optimal qx{s) giving the maximum value of Tc{H). If s is 
fixed in 1.0, we get the line with solid circles. 



We also study Tc{H) of an anisotropic spin-singlet 
state. As shown in our previous paperjj^ the linearized 
gap equation for the anisotropic spin-singlet state is 
given by 

A{x) = 2UT cos^ick^) ^ f da;'A(a;') 

xG"{x, x', ky, k^,u}n)G^x, x' , -ky, -fc^, -tj„). (13) 

where U is the nearest-site attractive interaction along 
the c axis. The energy gap is zero at the lines \ckz\ = 'k/2 
in this model. As in the isotropic case, eq. (13) is written 
as a matrix equation 



timal qxis), i.e., qxis) touches k^ ^ 
ky « ±7r/&. It may depend on the curvature of the Fermi 
surface and the Fermi velocity that which qx{s) gives the 
higher Hp'^^^ . We found that Hp'^^^ is the largest when 



kl,^ at fc, « or A^, = 2C/ ^ ^ ^ Df^, iky)K,^ (Q + NG)Ag , (14) 

I' N ky 

where Kky{qx) is given in eq. (12) and the coefficients 
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Dii,{ky) is defined by 



2M'^,U'{ky) 



(15) 



where ^ i'{ky ) is given by 

n 

X Jn^l+N{Vky.a)Jn-l'+N~j{Vky,a)- (16) 

In Fig. 3, we plot Tc{H) of an anisotropic spin- 
singlet superconductivity. Since the initial slope of the 
anisotropic spin singlet state is \/2 times larger than 
that in the isotropic stateP we take parameters as 
2ta ~ 2070K and 2tc ~ 4.0K in order to fit the ex- 
perimental results. Other parameters are same as in the 
isotropic pairing case. 
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Fig. 3. Transition temperature Tc{H) of the anisotropic pairing 
as a function of the magnetic field. 
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If [/ = A and tc ^ (r/^^^ = 0), T^iH) of an 
anisotropic spin-singlet is same as an isotropic pairing 
since S[^i,{ky) = Dj^i,{ky) = Sio6i'o6no- When tc ^ 0, the 
behavior of Tc{H) for an anisotropic spin-singlet state is 
different from that for the isotropic case, but the differ- 
ence between them is small. 

In conclusion, we have calculated Tc{H) of the 
isotropic and anisotropic spin-singlet superconductiv- 
ity in QID electrons. Although the Zeeman split- 
ting strongly suppresses the superconductivity, the crit- 
ical magnetic field is enhanced by choosing the opti- 
mal "nesting vector" and taking account of the higher 
harmonic terms in the energy dispersion. This re- 
sult seems toJie^nsistent with Tc{H) observed in the 
(TMTSF)2x|'BE3' which may show not only the pos- 
sibility of spin-singlet pairing in this system but impor- 
tance of higher harmonic terms. 
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